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Abstract 

The induced pseudotensor constant (weak electricity) of the nucleon is cal- 
culated in the framework of the chiral quark soliton model. This quantity 
originates from the G-parity violation and hence is proportional to m u — md- 

We obtain for m u — = —5 MeV a value of gr/gA = —0.0038. 
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I. INTRODUCTION 



The neutron /3-decay is a powerful tool to probe the structure of the nucleon. In partic- 
ular, it provides a precise measurement of the triplet axial constant of the nucleon g^, so 
that it is a touchstone for any model of the nucleon structure. The underestimation of the 
nucleon axial charge in the solitonic picture of the nucleon was for a long main critical point 
of soliton models of the nucleon. Recently it was shown [|T],|2|] that in the chiral quark-soliton 
model of the nucleon (xQSM) the rotational 1/N C corrections to the qa bring its value close 
to the experimental one. Also these corrections improve considerably the agreement of the 
electromagnetic characteristics of baryons ||-[5| calculated in the xQSM with an experiment. 

In the present paper we investigate the other than gA axial characteristic of the nucleon - 
induced pseudotensor constant (weak electricity) of the nucleon g T . The neutron-to-proton 
transition matrix element of the axial current = uj^'j^d can be written in terms of three 
form factors: 

( p (p')\ J f,\ N (p)) = U P (P){9A1^15 + M - M + 9P%lh}u n {p) , q = p -p, (1) 

where M p (M n ) is the proton (neutron) mass and we use the convention of Bjorken and Drell 
for Dirac matrices and spinors. The axial-vector gA and pseudoscalar gp[] constants were 
extensively analyzed theoretically and measured in experiments, while less is known about 
the pseudotensor constant gx- The pseudotensor current has the opposite G-parity to that 
of the axial vector current and hence is proportional to the parameter of isospin symmetry 
breaking. There are two different sources of isospin symmetry breaking: Electromagnetic 
interactions and u and d quark mass difference. In this work we calculate the hadronic part 
of the gT proportional to m u — uid in the limit of a large number of colors, N c — » oo. 

Even though in reality N c = 3, the limit of large N c furnishes a useful guideline. At 



l The pseudoscalar axial constant gp was calculated recently in ref. || in the framework of the 
chiral quark soliton model 
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large N c the nucleon is heavy and can be viewed as a classical soliton of the pion field 
An example of the dynamical realization of this idea is given by the Skyrme model 
A far more realistic effective chiral lagrangian of the xQSM is based on the interaction of 
dynamically massive constituent quarks with pseudo-Goldstone meson fields. It is given by 
the functional integral over the quark (ip) in the background pion field P|-|T2]]: 



exp (iS cS [n(x)]) = JvipVijexp^i J d 4 xi/;Di/?J , (2) 
where D is the Dirac operator 

D = ip-m-MW 5 . (3) 

U 15 denotes the pseudoscalar chiral field 

= expi7r a r a 7 5 = ^f^U + ^—^Ul (4) 

The m is the matrix of the current quark masses in = diag(m„,md). The M stands for 
the dynamical quark mass arising as a result of the spontaneous chiral symmetry breaking. 

The effective chiral action given by eq. (|2|) is known to contain automatically the Wess- 
Zumino term and the four-derivative Gasser-Leutwyler terms, with correct coefficients. 
Therefore, at least the first four terms of the gradient expansion of the effective chiral 
lagrangian are correctly reproduced by eq. @, and chiral symmetry arguments do not leave 
much room for further modifications. Eq. (||P has been derived from the instanton model of 



the QCD vacuum [fT2 |, which provides a natural mechanism of chiral symmetry breaking and 
enables one to express the dynamical mass M and the ultraviolet cutoff A intrinsic in eq. (0) 
through the Aq C d parameter. It should be mentioned that eq. (0) is of general nature: one 
can use eq. (0) without referring to the instantons. 

An immediate implication of the effective chiral theory eq. (0) is the quark-soliton model 
of baryons [T3J. According to these ideas the nucleon can be viewed as a bound state of 
N c (=3) valence quarks kept together by a hedgehog-like pion field whose energy coincides 
by definition with the aggregate energy of quarks from the negative Dirac sea. Such a 



semiclassical picture of the nucleoli is well justified in the limit N c — > oo - in line with more 
general arguments by Witten |7| . The further studies showed that the xQSM is successful in 
reproducing the static properties and form factors of the baryons using just one parameter 



set (see the recent review ||14|| ) . The powerful numerical method to carry out the calculation 



of the N — A splitting and nucleon matrix elements of arbitrary quark bilinear operators 



has been developed in refs. [J15| , |16|| . This method is also used in the present paper. 



II. COMPUTING WEAK ELECTRICITY 

The transition matrix element eq. ([J) can be computed as the Euclidean functional 
integral in the xQSM 

1 T T 

{PWl^d\N) = — hm^exp (ip - ~ Wo^) 



x J d 3 xd 3 yexp (-ip' ■ y + ip- x) J J Dip J Vtp^ 



(5) 



x J p (y, T/2) vrftfsd Jl(x, -T/2) exp \J d 4 z^D^ 
The nucleon current Jn (N = p, n) is built of N c quark fields: 

Jn(x) = J^ e ii---i N SjJ3TT3^aih( X )'''i J a Nc i Nc ( x )- (6) 

«i • • 4 0cpf c denote spin-flavor indices, while if • -in c designate color indices. The matrices 
^ ( jj 3 tt 3 c are taken to endow the corresponding current with the quantum numbers JJ 3 TT 3 . 

In the large N c limit the integral over Goldstone fields U in eq. (§|) can be calculated 
by the steepest descent method (semiclassical approximation). The corresponding saddle 
point equation admits a static soliton solution, an example of which is the hedgehog field 
configuration: 

U s (x) = exp [in ■ rP(r)]. (7) 

The P(r) denotes the profile function satisfying the boundary condition P(0) = tt and 
P(oo) = 0, which is determined by solving the saddle point equations (for details see 



Ref. [jTJJ]). The soliton is quantized by introducing collective coordinates corresponding 
to SU(2)j isospin rotations of the soliton (and simultaneously SU(2) spin in spin space): 

U(t,x) = R(t)U s (x)R\t), (8) 

where R(t) is a time-dependent SU(2) matrix. The quantum states arising from this quan- 
tization have the quantum numbers corresponding to the nucleon and A. 

Calculating the functional integral eq. (|5]) we obtain the following expression for the 
neutron to proton transition element of the axial current: 

(P|«7 M 7 5 d|iV> = N c (M p + M n ) J d 3 xe l ^J dR(j)* p {R) 

x /firf(x| — -— - 1 —— 7o7 , 75 r 1+a R\x)) <j> n (R), (9) 

J 2n \ to + iH + i(m u — m d )R^T A R J 

(~) (~) 

where 0|^(-R) is the rotational wave function of the nucleon (<f> p = 4>A, <p n = 0i 2 _i) given 
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by the Wigner finite-rotation matrix P,|T3 



47t 3 (R) = V2STl(-l) T+T ^ s _=l S3 (R), (10) 

and the integral over SU(2) group is normalized by / dR = 1. The one-particle Dirac 
Hamiltonian if in a background of the static pion field eq. (0) has a form 

H = 7 V(9 fc + tM^Uf + ^(m u + m d ) . (11) 

Projecting the general expression eq. (§) onto the pseudotensor structure one obtains: 

0L = Nc j^jLj dRm) 



x 



J 2tt \ u) + iH + i{m u — md)R'T' i R J 



Let us now show that the above expression is zero in the isospin symmetry limit (m u = 
md)- To prove this we introduce the following unitary transformation of the Dirac and Pauli 
matrices connecting them to the transposed ones: 

W^W- 1 = 7J, WW" 1 = -(r a ) T . (13) 
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Evidently then WHW 1 = H T . Using properties of the trace tr(M T ) = tr(M) and 
tiiWMW- 1 ) = tr(M) one can write: 

1 



tr (x| 



tr W(x 



~ tr ( x l 



ijj + iH + i{m u — m^Wr^R 
1 



75j R t r 1+i2 J R|x>^ = 



LO 



+ iH + i(m u — md)R^r 3 R 



75 



R^r^R^W- 1 



■ 75 R)r 1+i2 R\x) 



ijj + iH — i(m u — md)R^T z R 
This immediately implies that the pseudotensor constant given by eq. (|12j) is zero in the 
isosymmetric limit and first non-zero result appears expanding eq. (|T2"D in m u — rrid to linear 
order. The result for the pseudotensor constant in the leading order of 1/N C expansion 
(g? ~ iV c ) and the linear order in m u — has a form ]T7|] : 



9t 



iN c (m u - m d ) [ dw ^, 1 



■Tor 1 



M p + M n 24 J 2vr u + iH u + iH 

x e ab3 [ dR<p;(R)V^ a (R)V^(R)MR) ■ 



(14) 



The integral over soliton orientations in the second line of eq. ([14]) can be easily calculated 
by using the relations 



e ab3 V { ^ a (R)V^(R) = ±iV^ 3 (R) 



(1) 



and 



J dR (j)*(R)V^l i2 3 (R)(f) n (R) = -2/3. 



(15) 



(16) 



The functional trace in the first line of eq. (|T4l ) was estimated in ref. |17| by means of the 
gradient expansion: 

N c (m u -m d ). 



9t 
M p + M n 



9 ■ 96tt 2 M 



-Im d x 



-e klm e abm tr (r b d k U(r a U - tfV)$E7+) 



(17) 
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This approximation is justified only for a soliton of large size RM 3> 1. The real nucleon 
has a radius of order 1/M and hence the eq. (|H]) can be used only as an order of magnitude 
estimate. 



III. NUMERICAL RESULTS AND CONCLUSION 



In order to evaluate exactly the functional trace in eq. (|14T), we diagonalize the Hamil- 
tonian H eq. ( |iT|) numerically in the Kahana-Ripka discretized basis |L8 |. The constituent 



quark mass M is fixed to 420 MeV in our model by reproducing best many static baryon 
observables and form factors in the model (in particular, the isospin mass splittings for octet 
and decuplet baryons |19| , |T4|1 ). To make sure of the numerical calculation, we compare our 
results for g? with the analytical ones of the gradient expansion eq. (|17|) justified in the limit 
of large soliton size. Our numerical procedure is in good agreement within a few percent 
with the analytical results of the gradient expansion in the large soliton size limit. 

The results of our calculation are summarized in Table I. For completeness we give in 
Table I also results for gA obtained in ||. Let us note that the present result is comparable to 
a recent calculation of the nucleon pseudotensor constant with the QCD sum rule technique 
2(1 which gives gxj 9a = —0.0151 ± 0.0053. Both the QCD sum rule result and ours are in 



agreement with the bag model calculation (gr/gA — —0.00455) pO| , pi| , whereas they seem to 
be smaller than preliminary experimental data ranging from — 0.21±0.14 to 0.14±0.10. 
However, the accuracy of the experiment is not enough to be compared in a reasonable way 
with the results of theoretical calculations. 
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TABLES 

TABLE I. Axial vector qa and pseudotensor gj- constants of the nucleon as a function of the 
constituent quark mass M, m u — = —5 MeV. 

M [MeV] g A @ g T /g A 

370 1.26 -0.0029 

400 1.24 -0.0035 

420 1.21 -0.0038 

450 1.16 -0.0040 
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